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^ 1 ABSTRACT 

We calculate the 0{a ew m^ b ^/m'^ v ) and 0(a ew m^ b ^/m^ v ) supersymmetric elec- 



troweak corrections to the cross section for the charged Higgs boson production 
in association with a top quark at the Tevatron and the LHC These corrections 
arise from the quantum effects which are induced by potentially large Yukawa cou- 
plings from the Higgs sector and the chargino-top(bottom)-sbottom(stop) couplings, 
neutralino-top(bottom)-stop(sbottom) couplings and charged Higgs-stop-sbottom 
couplings. They can decrease or increase the cross section depending on tan/3 but 
are not very sensitive to the mass of the charged Higgs boson for high tan [3. At 
low tan/3(= 2) the corrections decrease the total cross sections significantly, which 
exceed —12% for m H ± below 300GeV at both the Tevatron and the LHC, but 
for m H ± > 300GeV the corrections can become very small at the LHC. For high 
tan/3(= 10,30) these corrections can decrease or increase the total cross sections, 
and the magnitude of the corrections are at most a few percent at both the Tevatron 
and the LHC. 

PACS number: 14.80.Bn, 14.80. Cp, 13.85. QK, 12.60. Jv 



1. Introduction 



There has been a great deal of interest in the charged Higgs bosons appearing in 
the two-Higgs-doublet models(THDM)[l], particularly the minimal supersymmetric 
standard model(MSSM)[2], which predicts the existence of three neutral and two 
charged Higgs bosons h,H,A, and H ± . When the Higgs boson of the Standard 
Model(SM) has a mass below 130-140 Gev and the h boson of the MSSM is in the 
decoupling limit (which means that H* 1 is too heavy anyway to be possibly produced), 
the lightest neutral Higgs boson may be difficult to distinguish from the neutral 
Higgs boson of the standard model(SM). But charged Higgs bosons carry a distinctive 
signature of the Higgs sector in the MSSM. Therefore, the search for charged Higgs 
bosons is very important for probing the Higgs sector of the MSSM and, therefore, 
will be one of the prime objectives of the CERN Large Hadron Collider (LHC). At 
the LHC the integrated luminosity is expected to reach L = lOOfb^ 1 per year in the 
second phase. Recently, several studies of charged Higgs boson production at hadron 
colliders have appeared in the literature[3,4,5]. For a relatively light charged Higgs 
boson, m H ± < m t — m b , the dominate production processes at the LHC are gg — > tt 
and qq — > tt followed by the decay sequence t — > bH + — > br + u T [Q]. For a heavier 
charged Higgs boson the dominate production process is gb — > tH~ [7,8,9]. Previous 
studies showed that the search for heavy charged Higgs bosons with m#± > m t + mj 
at a hadron collider is seriously complicated by QCD backgrounds due to processes 
such as gb — > ttb,gb — > ttb, and gg — > ttbb, as well as others process[8]. However, 
recent analyses [10, 11] indicate that the decay mode H + — > r + u provides an excellent 
signature for a heavy charged Higgs boson in searches at the LHC. The discovery 
region for H ± is far greater than had been thought for a large range of the (m#± , tan f3) 
parameter space, extending beyond m H ± ~ ITeV and down to at least tan/? ~ 3, 
and potentially to tan/5 ~ 1.5, assuming the latest results for the SM parameters 
and parton distribution functions as well as using kinematic selection techniques and 
the tau polarization analysis[ll]. Of course, it is just a theoretical analysis and no 
experimental simulation has been performed to make the statement very reliable so 
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far. 

The one-loop radiative corrections to H~t associated production have not been 
calculated, although this production process has been studied extensively at tree- 
level[7,8,9]. In this paper we present the calculations of the 0(a ew m^/m^ v ) su- 
persymmetric(SUSY) electroweak corrections to this associated H~t production pro- 
cess at both the Fermilab Tevatron and the LHC in the MSSM. These corrections 
arise from the quantum effects which are induced by potentially large Yukawa cou- 
plings from the Higgs sector and the chargino-top(bottom)-sbottom(stop) couplings, 
neutralino- top(bottom)-stop(sbottom) couplings and charged Higgs-stop-sbottom 
couplings which will contribute at the 0{a ew m^/m\ v ) to the self-energy of the 
charged Higgs boson. In order to get a reliable estimate this process has to be merged 
with the related gluon splitting contribution gg — > H~tb. This leads to a suppres- 
sion by about 50% at LO[12]. However, the complete one- loop QCD corrections are 
probably more important, but not yet available. 



2. Calculations 

The tree-level amplitude for gb — > tH~ is 

M Q = M^+M^\ (1) 

where and represent the amplitudes arising from diagrams in Fig. 1(a) and 
Fig. 1(6), respectively. Explicitly, 



M (s) = W9s -ufa) [2m t cot /3p£P L + 2m b tan /3rfP R - m t cot (3Y fcP L 

V2m w (s — m b ) 

-mbt&nP^fiPnjuip^e^T", (2) 

and 

M ® = 7k 5T "(Pt) t 2m * cot Prf p L + 2m b tan PPtPR - m cot /3>f JkP L 

V 2mw{t — mf) 

-m b twP>y' i / kP R ]u( Pb )e li (k)T? j , (3) 
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where T a are the SU (3) color matrices and s and i are the subprocess Mandelstam 
variables defined by 

s = (p b + k) 2 = (pt+PH-f, 

and 

t = (pt ~ kf = (p H - - p b f. 

Here the Cabbibo-Kobayashi-Maskawa matrix element VcKAi[bt] has been taken to 
be unity. 

The SUSY electroweak corrections of order 0(a ew m 2 , b Jm 2 v ) and 0(a ew mt b Jm^ v ) 
to the process gb — > H~t arise from the Feynman diagrams shown in Figs.l(c)-l(v) 
and Fig.2. We carried out the calculation in the t'Hooft-Feynman gauge and used 
dimensional reduction, which preserves supersymmetry, for regularization of the ul- 
traviolet divergences in the virtual loop corrections using the on-mass-shell renor- 
malization scheme [13], in which the fine-structure constant a ew and physical masses 
are chosen to be the renormalized parameters, and finite parts of the counterterms 
are fixed by the renormalization conditions. The coupling constant g is related to 
the input parameters e,m w , and mz by g 2 = e 2 /s 2 v and s 2 w = 1 — m 2 w /m 2 z . The 
parameter (3 in the MSSM we are considering must also be renormalized. Following 
the analysis of ref.[14], this renormalization constant was fixed by the requirement 
that the on-mass-shell H + lvi coupling remain the same form as in Eq.(2) of ref. [14] to 
all orders of perturbation theory. Taking into account the 0(a ew m 2 , b Jm 2 v ) Yukawa 
corrections, the renormalized amplitude for the process gb —> tH~ can be written as 

M ren = M W + M W + 5M Vl{s) + 5M Vl(t) + <5M s(s) + 5M s{t) + 5M V2{s) 

+5M V2(t) + 5M b{s) + 5M m = M,S s) + M (t) + ]T 5M l , (4) 

i 

where 5M V ^ S \ 6M Vl ®, 5M S ^ S \ 5M s{t \ SM V ^ S \ 8M V ^\ 5M b{ - s \ and 5M b ^ represent 
the corrections to the tree diagrams arising, respectively, from the gbb vertex dia- 
gram Fig.l(c)-l(d), the gtt vertex diagram Fig.l(f)-l(g), the bottom quark self-energy 
diagram Fig.l(i), the top quark self-energy diagram Fig.l(k), the btH~ vertex dia- 
grams Figs.l(m)-l(n) and Figs.l(p)-l(q), including their corresponding counterterms 
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Fig. 1(e), Fig. 1(h), Fig.l(j), Fig. 1(1), Fig.l(o), and Fig.l(r), and the box diagrams 
Figs.l(s) — 1(f). J2i $M l then represents the sum of the contributions to the Yukawa 
corrections from all the diagrams in Figs.l(c)-l(v). The explicit form of SM l can be 
expressed as 

+fi 2 P^P R }u(p b )e,(k), (5) 

where the C l are coefficients that depend on s,i, and the masses, and the f\ are 
form factors; both the coefficients C l and the form factors f\ are given explicitly in 
Appendix A. The corresponding amplitude squared is 

£|M ren | 2 = ]T|M (s) + M (t) | 2 + 2 J ReB(E W)(M« + (6) 

i 

where 

S|M W + M«| a = 1 2 , 2 [(m 2 cot 2 /3 + m 2 tan 2 /3)(p 6 -A;p t -A; 

^ 2N c m$ v (s-mfy 



m 2 b p t -k + 2p b - kp b ■ p t - m 2 b p b ■ p t ) + 2m\m 2 t {p b ■ k - m 2 )] 
1 

(i-m 2 ) 2 

- m 2 p b -p t ) + 2m 2 b m 2 (pfk-m 2 )] + 



+ — — 2 ^ [(m 2 cot 2 /3 + m 2 tan 2 /3) (p b ■ kp t ■ k + m 2 p b • fc 

1 



(s - m 2 ){t - ml) 
x [(m 2 cot 2 /3 + m\ tan 2 /3) (2p b • kp t ■ k + 2p b • ■ p t - 2(p b ■ p t ) 2 

- m 2 b p t ■ k + m 2 p b ■ k) + 2m 2 b m 2 t {p t ■ k - p b ■ k - 2p b ■ pt)}} , (7) 



V5M'(M (s) )t = --— J*f* C^h^fl (8) 

^ K ' Q4N C x lQn 2 m 2 w {s - m 2 b ) {r[ 1 Ul yJ 



and 



E^(^) f = - 7777 t/tS r^E^fl ( 9 ) 

^ 64A^c x 167r 2 m^(t - mf) i=1 

Here the color factor Nc = 3 and and hf^ are scalar functions whose explicit 
expressions are given in Appendix B. 
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The cross section for the process gb — > tH is 



167TS 2 



E\M ren \ 2 dt (10) 



with 



2 2 A I 

*min = ^ 2^ S ~ ^ + m ^-) )( s_ {m t -m H -) 2 ), 

and 

mf + — s 1 



tmax = g h-y(s- (m t + m^-) 2 )(s- (mt-m^-) 2 ). 

The total hadronic cross section for pp — > gift — > can be obtained by folding the 
subprocess cross section a with the parton luminosity: 

a(s) = C dz^a(gb^tH- at s = z 2 s). (11) 

Here ^/s and are the CM energies of the pp and gb states , respectively, and dL/dz 
is the parton luminosity, defined as 

~ d ~ = 2z J z2 ^fb/p(x,fi)fg/ P (z 2 /x,fi), (12) 

where fb/p(x, /i) and f g /p(z 2 /x, ji) are the bottom quark and gluon parton distribution 
functions. 



3. Numerical results and conclusion 

In the following we present some numerical results for charged Higgs boson pro- 
duction in association with a top quark at both the Tevatron and the LHC. In 
our numerical calculations the SM parameters were taken to be m w = 80.41C7el/, 
m z = 91.187GreV, m t = 176GeV, a s {m z ) = 0.119, and a ew {m z ) = j^-g[15]. And 
we used the running b quark mass ~ 3GeV and the one-loop relations[16] from the 
MSSM between the Higgs boson masses m h ^ Ht A,H± and the parameters a and (3, and 
chose m H ± and tan/5 as the two independent input parameters. And we used the 
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CTEQ5M[17] parton distributions throughout the calculations. Other MSSM param- 
eters were determined as follows: 

(i) For the parameters Mi,M 2 , and fi in the chargino and neutralino matrix, we 
put M 2 = 300GeV and then used the relation M 1 = (5/3)(#' 2 /g 2 )M 2 ~ 0.5M 2 [2] 
to determine M 1 . We also put /i = — lOOGeV except the numerical calculations as 
shown in Fig.6(b), where /i is a variable. 

(ii) For the parameters mj~ q ^ and A tib in squark mass matrices 



to simplify the calculation we assumed m| = m 2 ^ = m 2 ^ and A t = A b , and we put 
rriQ = 500GeV and A t = 200GeV. But in the numerical calculations of Fig. 6(a) 
A t = A b are the variables. 

Some typical numerical calculations of the tree-level total cross sections and the 
O^e^m^/m^,) SUSY electroweak corrections as the functions of the charged Higgs 
boson mass, A t = A b and /z, respectively, for three representative values of tan (5 are 
given in Figs. 3-6. 

Figures 3(a) and 4(a) show that the tree-level total cross sections as a function of 
the charged Higgs boson mass for three representative values of tan/3. For m^± = 
200Gel / the total cross sections at the Tevatron are at most only 0.7 fb and 0.1 fb 
for tan (3 = 2, 30 and 10, respectively, and for m H ± = 300GeV the total cross sections 
are smaller than 0.15 fb for all three values of tan/3. However, at the LHC the total 
cross sections are much larger: the order of thousands of fb for m H ± in the range 100 
to 240GeV and tan/3 = 2 and 30; and they are hundreds of fb for the intermediate 




(13) 



with 



m ll = m % + m l + m \ cos 2P(Iq L ~ e q sin 2 9 W ), 
M RR = rrijy p + m q + m z cos 2(3e q sin 6 W: 




(14) 
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value tan/3 = 10. When the charged Higgs boson mass becomes heavy(< 500 GeV), 
the total cross sections still are larger than 100 fb and 10 fb for tan/5 = 2,30 and 
10, respectively. For low tan/3 the top quark contribution is enhanced while for high 
tan (3 the bottom quark contribution becomes large. These results are smaller than 
ones given in ref.[8,9] because we used the running b quark 3GeV in the 

numerical calculations. We have confirmed that if we chose nib = 4.5GeV, our results 
will agree with ref.[8,9]. 

In Figs. 3(b) and 4(b) we show the corrections to the total cross sections relative 
to the tree-level values as a function of m H ± for tan/3 = 2, 10, and 30. For tan/3 = 2 
the corrections decrease the total cross sections significantly, which exceed —13% 
for m H ± below 300GeU at the both Tevatron and the LHC. But the corrections 
decrease as m#± increase. For example, as shown in Fig. 4(b), the corrections range 
between —13% ~ 0% when m H ± increase from 300GeV to ITeV at the LHC. For 
high tan/3(= 10, 30) these corrections become smaller, which can decrease or increase 
the total cross sections depending on tan /3, and the magnitude of the corrections are 
at most a few percent for a wide range of the charged Higgs boson mass at both the 
Tevatron and the LHC. 

In Fig. 5 we present the Yukawa correction from the Higgs sector and the genuine 
SUSY electroweak correction from the couplings involving the genuine SUSY par- 
ticles(the chargino, neutralino and squark) for tan (3 = 30 at the LHC, respectively. 
One can see that the Yukawa correction and the genuine SUSY electroweak correction 
have opposite signs, and thus cancel to some extent. The former decrease the total 
cross sections, which can range between —8% ~ —4% for m H ± below 300GeU, but 
the latter increase the total cross sections, which range between 10% ~ 7% for m H ± 
in the same range. In such a case the combined effects just are about 2% ~ 3%. 

Figs. 6(a) and 6(b) give the corrections as the functions of A t = A b and // for 
m H ± = 300GeU at the LHC, respectively, assuming tan/3 = 2,10 and 30. From 
Figs. 6(a) and 6(b) one sees that the corrections increase or decrease slowly with 
increasing A t = A b and the magnitude of //, respectively, for tan (3 = 30, 10, and the 
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corrections are not very sensitive to both A t = Ab and \i for tan/5 = 2, where the 
corrections are always about —12% and —13%, respectively. In general for large values 
of A t and small values of tan (3 or large values of \i and tan /3, one can get much larger 
corrections since the charged Higgs boson-stop-sbottom couplings become stronger. 
For tan (3 = 30, comparing Fig. 4(b) with Fig. 6(b), we can see that the corrections 
indeed become larger as the values of \x increase. But for tan/3 = 2 from Fig. 4(a) and 
Fig. 6 (a) we found that the corrections almost have no change when A t = A b become 
larger. Obviously the effects from the stronger couplings have been suppressed by the 
decoupling effects because with an increase of A t = A b all the squark masses are still 
heavy, which almost is same as discussed in Ref.[18]. 

In conclusion, we have calculated the 0(a ew m^/m^ v ) and 0(a ew m^/m^ v ) SUSY 
electroweak corrections to the cross section for the charged Higgs boson production in 
association with a top quark at the Tevatron and the LHC. These corrections decrease 
or increase the cross section depending on tan f3 but are not very sensitive to the mass 
of the charged Higgs boson for high tan (3. At low tan/3(= 2) the corrections decrease 
the total cross sections significantly, which exceed —12% for m H ± below 300GeV at 
both the Tevatron and the LHC, but for m H ± > 300GeU the corrections can become 
very small at the LHC. For high tan/?(= 10,30) these corrections can decrease or 
increase the total cross sections, and the magnitude of the corrections are at most a 
few percent at both the Tevatron and the LHC. 
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China, the Doctoral Program Foundation of Higher Education of China, the Post 
Doctoral Foundation of China, a grant from the State Commission of Science and 
Technology of China, and the U.S. Department of Energy, Division of High Energy 
Physics, under Grant No.DE-FG02-91-ER4086. S.H. Zhu also gratefully acknowledges 
the support of the K.C. Wong Education Foundation of Hong Kong. 
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Appendix A 



The coefficients C l and form factors f\ are the following: 

2 2 2 

s) = m b (jVxit) _ m t qs(s) nr ~ 



rr$y(s-miy m 2 w {i - m 2 )' m 2 w (s - m\) 2 ' 



,2 



t) = m t c V 2 (s) = m b m t c v 2 (t) = m b m t 

m 2 w (i-m 2 ) 2 ' ' m 2 w {s-m 2 h Y ' m 2 w (i - m 2 ) ' 

s) = C b(t) = m tm b 



rrv 



w 



= V {1) [M9^ S) -9^ S) )-2 Pb -k g ^% 

= V {2) [rn b (g^ s) -g^ s) )-2 Pb .kg^% 

= T,V[2(g™ + ^ 1(S) ) + ™ b (g^ + ^ 1(S) ) + 2p b • rf (s) ], 

= t/W^ + g™) + m b (g^ + g? (s) ) + 2 Pb ■ kg^% 

= v^-(g^ s) + g^ s) )+Mg^ s) + g^ s) )i 

= ^[-(^ (s) + ^ (s) ) + ^(^ (s) + ^ (s) )], 

= v { %: i{s) + 2 g ^ + rn b (g^ s) -g^ s) )l 

= vV s) + 2g? {s) + Mg? {s) -g? { n, 

= 2(?r 2(s) + 2m t cot PiSA? + 8Af + SA^) - 2m t g^ (s) + 2m b g? {s \ 

= 2g^ s) + 2m b t a nf3(5A% + 6A$ + 5Ag } ) + 2mtf? w - 2m t g\^\ 

_ 1 f V 2 (s) rV 2 (s) _ 1 ,V 2 (s) 

= 2«?r 2( * ) + 2m, cot /3(5A« + 5Af + Mf) - 2m t gf {t) + 2m^ 4 y2(<) , 

= 2(7 2 y2W + 2m fe tan/3(5Ag ) + 5A^ + <5Ag } ) + 2m b g^ {t) - 2m t g? {t) , 



1 rV 2 (t) rV 2 (t) 1 

2 ■/ 5 ) is 2 



2^-M^r+^(^+^ro], 

27 7 ( 2 )p b -M# ) + m 6 (^ ) + ^ (s) )], 
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/ 3 S(S) = V 2 >M<7; W + g s 5 (s) ) + 2{m\ + Pb ■ k)gi s) + (ml + 2p b ■ k)g s 3 (s) + m 2 b gf s \ 
= 2^{m b (gf s) + gf s) ) + 2{m\ + p b • k)g s 2 {s) + m 2 b g s 3 (s) + {m\ + ^ • k)g?\ 

Ms) _ 1 ,s(s) f s(s) _ 1 ,s(s) 

/f W = E ? / f ) {?/ 2 )[2m fc (-3A3i2 + (l-Ci)A 7 ) + m?(A> + ^12-^22 

i=H°,h°,G°,A° 

-D 32 ) - m 2 m b (D 23 + 2D 39 ) - 2m bPb ■ k(2D 36 + D 24 + d(D + D 12 )) 
+2m bPt ■ k(D 25 + D 310 ) + 2m bPb ■ p t (D 26 + 2D 38 )] + rj^[2m t (-3D 313 + (1 
+(i)D 27 ) - m 3 t (D 33 + (1 + (i)D 23 ) + m 2 m t (D 13 - 2D 38 + (1 + &)(A) 
-L> 22 )) + 2m t p fe • k(D 13 - D 3W - (1 + Q{D 12 + L> 24 )) + 2m tPi • k(2D 37 
+ (1 + 0)D 25 ) + 2m tPb ■ p t {2D 39 + (1 + 0)^ 26 )]} 
(-k, -p b , p t , m b , m b , m u m t ) 

-^tt E NkAN^RityRjtyaijDzri-k, -p b ,Pt, m^m^m^m^), 

P i,j,k 

f Ks) = tf%W„ v W,L l ~R l ,N kl ~Nti, 

fs is) = E V? ) W 1) [-4D 27 + 2ml(D 22 -D -(l-( t )(D 12 + D 22 )) 

i=H°,h°,G°,A° 

+2m 2 (D 23 - (1 + Ci)D 26 ) + Ap t ■ k(D 26 - D 25 )} + r]^2m t m b (l + Q(D 22 

-D 12 - D 26 )}(-k, -p b ,p t ,m b ,m b ,mi,m t ) 

^ aij[—m t N ki Nl 3 Ri(b)Rj(t)D 2& + m b N* k4 N k3 L i (b)L j (t)(D 12 
sin ■ ■ 7 

+D 22 ) + m^oN^N^RiifyLji^Du^-k, -p b ,p t ,m b .,m b .,m^m~ t .), 
ff s) = E V? ) W 1) [l2D 313 + 2ml(2D 38 -D 13 + (l-C l )(D 13 

i=H°,h°,G°,A° 

+D 26 )) + 2m 2 t (D 33 + (1 + Q)D 23 ) + Ap b ■ k(D 25 + D 310 ) - Ap t ■ k(D 23 

+2D 37 ) - 4p t -p b (D 23 + 2D 39 )] + V i2 hm t m b {l + Q)(D 13 + D 23 

-D 26 )}(-k, -p b ,p t ,m b ,m b ,mi,m t ) 
8\^2mw 



E <Jij [-m t Nk4Nl 3 Ri(b) Rj (t) D 23 + m fe iV* 4 7V fc3 L i (6)L i (t)(L» 1 3 



+D 26 ) + m^ k N* kA Nl 3 Ri (b)Lj (t)D 13 \ (-k, -p b , p t , m~ b . , m~ b . , m^o , m~ t . ) , 
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/r (s) = E r ] l 3) {r ] {1) [HD 27 -D 311 )+mt(D 11 -2D 12 -2D 22 

i=H°,h°,G°,A° 

-2D 36 + (1 + Ci)(D + D 12 )) - m 2 (2D 23 + 2D 37 + (1 + Q)D 13 ) - 2p b ■ k(D 12 
+2D 24 + 2D 3A ) + 2p t ■ k(D 13 + 2D 25 + 2D 35 ) + 2p t ■ p b (D 13 + 2D 26 
+-D310)] + il {2) m t m h (l + Ci)(D 12 - D 13 - D )}(-k, -p b ,p t ,m b ,m b ,mi,m t ), 
f b s (s) = f b 7 (s \v (1) ~V { \ 

f b 9 (s) = E ^ 3) {^)2m,[- J D 13 - J D 26 + (l + (i)(D 12 + D 24 )} 

i=H ,h°,G°,A 

+7] {2) 2m h [-D 22 + L>24 + C*(A> + %D 12 + D 24 )]}(-k, -p b) p t ,m b) m b) m h m t ) 
_8vW ^ aijNk4N;3mRj{t){Di2 

bm Z ' J i,j,k 
+D 24 )(-k,-p b ,p t ,m~ b .,m~ b .,m i o,m i .), 

/if = E V? ) W 1) 2m t [D 23 -(l + Ci)D 25 ] -r]^2m b [-D 26 + D 25 

i=H0,h°,G°,A» 

+d(D 13 + £> 2 5)]}(-&, -p b ,p t ,m b ,m b ,mi,m t ) 
+D 25 )(-k,-p b ,p t , m b . , m b . , m^o , m t ~. ) , 

/i2 S) = yS -) fa (1 W 2 U~i*, 

where D , Dij, D^^ are the four-point Feynman integrals [19]. The explicit forms of 
<5M yi W,6M s W,6M 6 W can be respectively obtained from <5M Vl W, <5M S W, <5M fe M by 
the transformation U defined as 

v^vu s^i k^-k, d 3) -d 4) , i^-uf , 

m t ^m b , r] {1) ^ r] (2) , \ b ^\ t , mi.^m b ., U i2 ^ V* 2 , N i3 <-> iV* 4 , 
Li(6) <-> L,(t), i2i(6) <-> i^t), p£P L(R) <-> rf^(L), 7 M <-> 7 M 
All other form factors not listed above vanish. In the above expressions we 



have used the following definitions: 

R \. 

\f2m-w cos (3 ' v^mvi/ sin (3 



— m b tan (3, r/^ = m t cot (3, \ b = —= , A t — 
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cos6 q , L 2 (q) = -sm9 q , Ri{q) = sin6 q , Ri^q) = cos8 q , 



cos 2 a m sin 2 a 



„(i) _ bm a „(i) _ tan 2 _(i) _ i 

cos 2 p Vh ° ~ cos 2 p Va ° ~ ^' Vg ° ~ 

sin^a (2) _cos^ (2) _ 2 (2) _ 

sin 2 /?' V "sin 2 /T ^o-cot^, ^ - 1, 

(3 ) _ sin a cos a (3) _ (3) _ 

" sin /? cos F VG ° ~ ~ VA ° ~ ' 

: 2 /^J £g _ = 2' £ff- = 2 ^ al1 ^' = 2> 

tan 2 /5, ^ = 1, 4^=cot 2 /3, = 1, 

Cff° = Cft° = Cff- = ~(a° = ~(g° = ~(g- = 1) 

^ (sin 2/3 - ml tan g ± m * COt g (t) 
v 2 "^vk 

+-^(tan/? + cot /3)Ri(b)Rj(t) — (>Lt - i4 6 tan/3)i2 i (6)L J -(*) 

V2^iy v2mvi/ 

'"' (//- A t co\,p)Li{b)Rj(t), 



V2m 



w 

^ - 2C 24 + m 2 (-2C n + C 12 - C 21 + C 23 ) - s(C 12 

i=HO,h°,G°,A° 1 

+C23) ](-£>&, -k,mi,m b ,m b ) + [-F + F 1 + 2m 2 Gi 
-(1 + Ci)2m 2 Go](m 2 ,mi,m 6 )}, 

2 2{^ (1) [i - 2C 24 + m 2 C + m 2 (-C - 2C n + C 12 - C 21 + C 23 ) 

i=H-,G~ 1 

-s(Ci2 + C 23 )](-p 6 , -fc,m i ,mt,m t ) + [^(-^o + ^i) - 2m 2 0G 
+ C (3) )(Gi - 0G )](m 2 ,m 4 ,m,)} 



4m 2 



]T{A 2 [i? 2 (6)|iV, 3 | 2 (-Fo + Fx) + m 2 b \N i3 \ 2 (-G + d) - 2m b m^ 



xL,(b)R,(b)N* 3 2 G ](mlm h ,m^) + [^m^+X^L^R^VgU^Gv 
+A 2 i? 2 (t)|^ 2 | 2 (-F + Fx) + m 2 (A 2 it! 2 (0|^2| 2 + A 2 L 2 (t)|f/ l2 | 2 )(-G 



+G'i)](m 6 ,m t -.,m-_+) - 2A 6 -R J -(6)|7V" i3 | C 24 (-p 6 , -A;, m^o, mg., m~ b .) 
\ 2 R 2 (t)\V l2 \ 2 ( 

92 » &\v a - U* 2 , N t3 <-> TV* , L,-(6) - A 6 L,(t) <- \ t Rj(t)), 



-2A 2 J R 2 (t)|^ 2 | 2 C , 24 (-^,-A;,m-+,m, > m f3 )}, 



^ ^f^m^Co + 2Cn + C 2 i + Ci(C + Cn)](-p 6 , -k,mi,m b ,m b ) 

i=H°,h ,G°,A 
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2 

+ 4:m b [tf\c + 2Cn + C 21 ) + ^( i (C + Cn)](-Pb,-k,m l ,m t ,m t ) 

i=H-,G~ m b 

+^E{^K°^( fo )^( & )^ 2 (^o + C'n)-m f) L J 2 (6)|iV J 3| 2 (C' 11 
m b ij 

+C21)] (~Pb, -k, m^o , m~ b . , m~ b . ) 

+ [m-+A 6 A t L J (t) J R J (t)V;*2^(Co + C n ) - m 6 A^(f)|E/ i2 | 2 (C'ii 
+C 2 i)](-p b ,-k,m^+,m i .,m i .)}, 

^ l(s) (d 1} - d 3) , - ££, iV, 3 - iV£, Lj(b) <-> R 3 (b), \ b L 3 (t) <- \ t R 3 (t)), 
X ^^m^Co + Cn + (iC )(-p b , -k,mi,m b ,m b ) 

i=H°,h°,G a ,A° 

- X 2m t>[d 3) (Co + Cn) H |CiCo](-Pb, -k,mi,m t ,m t ), 

i=H-,G~ m b 

X 2 Vi 1] (Ci2 + C 23 )(-Pb, -k,mi,m b ,m b ) 

i=H ,h°,G°,A 

+ X 4^ 1) (C 12 + C 2 4){-Pb,-k,m i ,m t ,m t ) 

i=H-,G~ 

^{A^(6)|JV,3| 2 (C 12 + C 23 )(-ft ) 4 ) m^m^ m.) 



+ A 2 i? 2 (t) \V l2 \ 2 (C12 + C 23 )(-p 6 , — fc, m -+ , m t ~ . , m~ t . ) } , 
ff8 yi(s) (d 1} <- V2 <- E? 2 , iV, 3 <- iV* 3 , L,-(6) <- it^fe), \ b Lj(t) <- XtRjit)), 
X ^ (3) {^ (1) ["^ + 4^24 + m 2 t (C + 2C U + C*(Co + Cn) 

i=H°,h°,G°,A° 1 

+C 21 - C12 - C 23 ) + m|_ (C 22 - C 23 ) + s(Ci2 + C 23 )] + ?7 (2) Tn b rrit [QCn 

+ (1 + C»)Cb]}(-Pt, -pH-,mi,m t ,m b ) + 4 ^" 7 ^ J2[m t Ri(b)Rj(t)NZ 3 N k4 

sin - - 7 

x(-Cn + C12) + m^oLji^RiifyN^N^Co^iji-pt, -p H - ,m ± ^m~ ho mi^ 

J2 vPW^miCo + C n + d(C + C 12 )] + J^Ci™^} 
(-p t , -p H -,mi,m t ,m b ) 

— a ir,9/T ^ R i(b)Rj(t)N; 3 N k4 a i:j C 12 (-p t: -pH-.m^m^mi.), 

' i,j,k 

93 2(s \r} {1) <-> iJ (2 Ui ~ Ri,N kl ~ TV*), 
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E + 4C 24 + m 2 (C + 2C U + Ci(C + C u ] 

i=H ,h°,G°,A 

+C21 — C12 — C23) + (C22 — C23) + i(Cl2 + C23)] 
+77 (2) TO 6 m t [C + C«(Co + Cii)]}(-P&,Pjt-, "i 6 , m t ) 

' ink 



SU1 - 



+ m X° L i (*) ^ ( 6 ) N k3 N k4 C o\ <7ij ( -P6 , PH- , ™x° ' m bi ' m tj ) > 

E ^{^WCo + Cn + O(O) + C 12 )] + ? 7 (2) C^ t C' 12 } 

i=H°,h ,G ,A° 

(-p b ,PH-,mi,m b ,m t ) 

1 ink 



E m^f^-CiFo^ + fc, mi, m 6 ) + [O^o - 2m 2 (l + Ci)G 

o 

+2m 2 b G 1 }(ml,m i ,m b )} + E 2m 4 \d F o(Pb + k, rrii, m t ) 



i=H-,G~ m b 



2 

+ [-2m 2 0G + m 2 ^ + d 3) )(Gi - OGo) + cA F o]K, m 4 , m,)} 



ml 



Ei-^MmmtiFoiPb + k, m bj , m,o) 
-m-+A 6 A t L J (fe)i? J (fe)^C/*F (p 6 + k,m- t .,m^) + [m^|iV i3 | 2 (-Go + G h 
-m^X 2 b L j (b)R j (b)N;i(2m 2 b G - F )](m 2 , m 6j , m s? ) + K(A 2 L 2 (t)|f/ l2 | 2 
+A 2 i? 2 (t)|\/ i2 | 2 )(-G + Gi) - m^X^L^R^V^U^mlGo 
-F )}(m 2 b ,m { .,m^+)}, 

E ^(--^o + F l )(p b + k,m,i,m b ), 

i=H°,h°,G°,A° 

E ^Vo - F 1 - 2m 2 b G 1 + 2(1 + ( t )m 2 b G ](m 2 b ,m t ,m b ) 

i=H°,h°,G°,A° 

+ E 2{C (1) (-^o + i 71 i)(P6 + A;,m J ,m t )-[d 1) (-^o + i 71 i) 

i=H-,G~ 

-20m 2 G + m 2 fe (1) +d 3) )(Gi - CG )](m 2 , m„ m t )} 
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£{A 2 [i? 2 (&)|Ar j3 | 2 (-F + Fx) + \N t3 \ 2 m 2 b (-Go + G 1 ) 
-2m h m%Lj (b)Rj (b)N* 2 G ] {m 2 b , m~ b . , m^o) 

+ [XlR](t)\V l2 \\-F Q + F ± ) + m 2 b (\ 2 R 2 (t)\V t2 \ 2 + A^ 2 (t)|C/ i2 | 2 )(G' 1 - G ) 
-2m b m~+L j (t)R j (t)X b XtV^U* 2 G ](ml, m f ., m~+) 
-A^(&)|JV i3 | 2 (-*b + F!)(p 6 + fc,m^.,m^) 
- A 2 R) (t) \V t2 \ 2 (-F + F 1 )(p b + k, m h , m -+ ) } , 

^3 (s) (d 1} <-> d 3) , Va - 0£, ^ 3 - iV£, L,(6) <-> R,(b), \ b L 3 (t) <-> \ t R 3 (t)), 

g *\N£ _> ^3, - vb, c/* - c/ i2 ), 

LY "-(1 - cot 2 9 W ) [2m 2 (In ^ - 1) + m 2 + m 2 - ^m 2 w + m 2 b F 



3m 2 / /l tv /i 2 7 1 6 



+ (m 2 - m 2 - 2m 2 w )F 1 ](m 2 w ,m b ,m t ) + cot 2 e w {--[(g b v ) 2 + (^) 2 

om w o 

+ (^) 2 + (<&)>! + Wv) 2 + (^) 2 ) (2m 2 In -± + m 2 F - 2m|F 1 ) 

o 

-((^) 2 - (^) 2 )3m 2 F ](m|,m t ,m t ) + [{{g b v f + («£) 2 )(2m 2 ln^ 

/i 

+m 2 F - 2m|Fx) - {{g h v f - (g b A f)3m 2 b F ](m 2 z ,m b ,m b )} + -^[(cot 2 /? 
-l)m 2 F + (m 2 - m 2 b - 2m 2 cot 2 /3)Fi + (m 2 cot 2 /3 + m 2 tan 2 /3 
+2m 2 )m 2 G — (m 2 cot 2 /3 + m 2 tan 2 (3)m 2 H - G±] (m 2 H - , m t , m b ) 
+ E -F - 2m 2 (G + 0G - G 1 )](mlm i ,m b ) 

i=H°,h°,G°,A° w 



2„(2) I 



-m t ry- 



-(1 + 2Ci)F + Fi + 2m 2 (l + 0)G - 2m 2 G 1 ](m 2 ,m i ,m t )} 

1 ' ,2r t (l)/ jp_ , Z7 1 \ o^2 A ^_ , 2/ t (l) , t (3)> 



+ E ^W 2 [d 1) (-^o + Fi)-2m 2 C i Go + m 2 (d 1) +d 3) ) 



i=H-,G~ m W 



x(G 1 - C^oMm 2 ,™^) " m 2 [-^C^o + efV^o + F) + 2m 2 0G 
-™ 2 (d 2) + d 4) )(Ci - CiG )]K, mi, m 6 )} - 2N c ^{^ijG 

1 2 2 

+Ri(b)Rj(t)m t m b (taii 2 f3 — cot 2 /3)]}(m 2 ^_, m t - , ^.), 

-2 E{\ 2 ["— ^(0^-(*)^4 (^o - ^ 2 Go) + \N ti \ 2 (R 2 (t)(-Fo + Fx) 
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2m~+ 

-m 2 (-G + G 1 ))](m 2 t ,m- tj ,m jc o) + [ ^\ b \ t L 3 {b)R 3 {b)U* 2 V*(F 

Tflf 

-m 2 t G ) + \ 2 b R 2 (b)\U t2 \ 2 (-F + Fx) - m 2 t (\ 2 L 2 (b)\V l2 \ 2 + \ 2 b R 2 (b)\U l2 \ 2 ) 
x (-G + Gi)](m t 2 , m b .,m-+)}, 
5A { ? = 2^{A 6 2 [|M3| 2 ( J R|(6)(-F + F x ) + m 2 b (-G + G,)) - 2m b m^L 3 (b)R 3 (b) 



xN* 2 G ](m 2 b ,m b m^o) + [-2m-+m b \ b \ t L 3 {t)R 3 {t)U* 2 V*G 



Xi 

+\ 2 b L 2 (b)\U l2 \ 2 (-F + F 1 )+m 2 b (\ 2 R 2 (m 2 \ 2 ^ 
x(-G + Gi)](m 2 ,m fj ,m-+)}, 
5A% = £ -V{m 2 r ? f ) [-F + F 1 -2m 2 (G + OGo-G 1 )](m 2 ,m J ,m,) 

-m 2 ^ (1) [-F + Fx - 2CiF + 2m 2 (l + Q)G - 2m 2 G 1 ](m 2 ,m i ,m b )} 
+ E ^K[d 2) (-^o + Fi)-2m 2 C i G + m 2 (C (2) +d 4) )(Gi 

i=H-,G- m W 

-QG^imlm^m,) - m 2 [-M Cij p + ^)(-F + F x ) + 2m 2 QG 

m b 

~m 2 M ] +C (3) )(Gi -QG )}(m 2 b , mi ,m t )}, 
*Ag> = 5Kf{U\ 5kf = 5kf(U). 

Here C , C^- are the three-point Feynman integrals[19] and C 24 = — + C 24 , while 

j?, \ f 1 i n, r -g 2 y(l - y) + m 2 (l -y)+ m 2 y 

F n {q,m 1 ,m 2 ) = / dyy In , 

J o /i z 

G (g mx m 2 ) = - C dy V n+ \ l ~ v) 

Jo -q 2 y(l - y) + m\{l - y) + m\y 

and 

t 1 4 • 2n t 1 6 1 2 • 2^ 6 1 

#v = 2 ~ 3 sm 9 A = 2' #v = - g + 3 sm °w, 9 A = - ^ 

which are the SM couplings of the top and bottom quarks to the Z boson. The 
definitions of 9 q , U i3 , V i3 , N i3 , /i, A q can be found in ref.[2]. 
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Appendix B 



12 



4m t ri ( - 2 \2pb ■ k — p® ■ p b ) — Am b ri < ^\p^ ■ Pt+Pt • k), 

2r] {2) {2p b ■ kp b ■ p t - m\p t ■ k - 2p {%) ■ p b p b ■ Pt) + 2m b m t r)^\p b ■ k 
-2pW-p b ), 

2r] (2 \m 2 p b ■ k - 2p {l) ■ p t p b ■ p t ) + 2m b m t m t 7] {1 \p t ■ k - 2p w • p t ), 

4^(2) ( p W . PbPt . k - p^ ■ kp b -pt-Ph- kp {i) ■ p t - 2p b ■ kp t ■ k) 
—A.m b mtrp~ 'p^ • k, 

Am t r^ > p b ■ k(pb ■ k — p^ ■ p b ) — Am b r]^p^ ■ p b pt • k, 
4% (1) ~77 (2) ), 

Am t r] {2) p b -k(pfk- p {l) ■ p t ) - Am b r] {1) p t ■ kp {t) ■ p t , 



where the index i represents the two channels s and t, and p^ = p b , p® 



Pt- 
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Fi gure 1: Feynman diagrams contributing to 0(a ew rn^^/m^y) Yukawa corrections to 
gb — > tH~: (a) and (6) are tree level diagrams; (c) — (v) are one-loop diagrams. The 
dashed lines represent H,h, A, ,G° and G ± for diagrams (c) and (/); H,h,A and G° 
for diagrams (m),(p),(t) and (u); t,b,H,h,A,H ± ,G° and for (i) and (fe), where the 
solid lines represent charginos and neutralinos if the dashed lines represent squarks. For 
diagrams (d) and (</), the solid lines in the loop represent x° and x + an d the dashed lines 
represent squarks. 




Figure 2: Self-energy Feynman diagrams contributing to renormalization constants: The 
dashed lines represent t, b, H, h, A, H^, G° and G ± for diagram (a), where the solid lines 
represent charginos and neutralinos if the dashed lines represent squarks. 
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Figure 3: The tree-level total cross sections (a) and relative one-loop corrections (b) 
versus m H ± at the Tevatron with y/s = 2 TeV. The solid, dashed and dotted lines 
correspond to tan/5 = 2, 10 and 30, respectively. 




Figure 4: The tree-level total cross sections (a) and relative one-loop corrections 
(b) versus m H ± at the LHC with y/s = 14 TeV. The solid, dashed and dotted lines 
correspond to tan/5 = 2, 10 and 30, respectively. 
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Figure 5: The radiative correction from top, bottom quarks (dashed line) and genuine 
SUSY particles (dotted line), as well as total contributions (solid line) when tan/3 = 
30 at the LHC with yfe = 14 TeV. 
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Figure 6: Relative one- loop corrections versus A t , Aj, (a) as well as \x (b) at the LHC 
with y/s = 14 TeV, where m H ± = 300GeV and the solid, dashed and dotted lines 
correspond to tan/3 = 2, 10 and 30, respectively. For (a), /i = — lOOGeV, and for (b), 
A t = A b = 200GeV. 



